The quantum brachistochrone problem for non-Hermitian Hamiltonians 
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Recently Bender, Brody, Jones and Meister found that in the quantum brachistochrone problem 
the passage time needed for the evolution of certain initial states into specified final states can be 
made arbitrarily small, when the time-evolution operator is taken to be non-Hermitian but VT- 
symmetric. Here we demonstrate that such phenomena can also be obtained for non-Hermitian 
Hamiltonians for which PT-symmetry is completely broken, i.e. dissipative systems. We observe 
that the effect of a tunable passage time can be achieved by projecting between orthogonal eigen- 
states by means of a time-evolution operator associated to a non-Hermitian Hamiltonian. It is not 
essential that this Hamiltonian is PT-symmetric. 



■ 

o ■ 
o : 

> ■ 

O: 
-xt- : 

(N ' 

> : 

in ■ 
(N ■ 

cn : 
o . 
r» ■ 
o ; 

^ ■ 

a,: 

^— > ■ 
c ■ 
a • 

=5 

cr 



X 



PACS numbers: 03.65.Xp, 03.65.Ca, 03.65.-w, 11.30Er 



I. INTRODUCTION 

To find the brachistochrone is one of the oldest prob- 
lems in classical mechanics tracing back to Newton and 
Leibniz. It consists of finding the trajectory between two 
locations of a particle, subject to a gravitational field, 
for which the transition time becomes minimal. This 
problem can be generalized to a relativistic [l[ and to a 
quantum mechanical setting 0, [H, 0, H @ • In the latter 
case one seeks the minimal time t =: r, referred to as 
passage time, such that 



-irH 



(1) 



for given initial and final states and \ipf), respec- 
tively. Equality can be achieved by possibly tuning some 
parameters in the Hamiltonian H . 

Bender, Brody, Jones and Meister [|| extended this 
treatment by allowing also non-Hermitian Hamiltonians 
in ([1]). The surprising result found in Q was that when 
involving non-Hermitian, but PT-symmetric Hamiltoni- 
ans, in the evolution operator, the passage time can be 
made arbitrarily small by varying a parameter in H while 
keeping the transition frequency between two states con- 
stant. At present this phenomenon is an observation and 
no explanation has been provided as to where this effect 
might originate from. 

One might suspect that one could make the PT- 
symmetry responsible for the observation and seek for 
similar arguments as those which allow to explain the re- 
ality of the spectrum of a 'PT-symmetric non-Hermitian 
Hamiltonian. See for instance [7|, |8|, l£j, [l0| for recent re- 
sults and reviews. The main purpose of this paper is 
to investigate whether PT-symmetry can be utilized as 
well to explain the observed phenomenon of a tunable 
passage time. In fact, we find that the same conclusion 
can be drawn when considering non-Hermitian Hamil- 
tonians with complex eigenvalues describing dissipative 



systems, i.e. those for which PT-symmetry is definitely 
broken. This means the possibility of arbitrarily small 
passage times results from the non-Hermitian nature of 
the Hamiltonian involved in the time-evolution operator 
and not its PT-invariance. 

Our manuscript is organised as follows: In section II we 
derive passage times for PT-symmetric Hamiltonians for 
various different types of initial and final states. In sec- 
tion III we perform a similar analysis for non-Hermitian 
Hamiltonians with complex eigenvalues deriving similar 
phenomena as in section II. We state our conclusions in 
section IV. 



II. PSEUDO HERMITIAN HAMILTONIANS 

We start by considering PT-symmetric or rather 
pseudo Hermitian Hamiltonians. We recall [TJ E2, EH 
QUI that a non-Hermitian Hamiltonian H is said to be 
a pseudo-Hermitian operator, if there exists a Hermitian 
operator ry, with regard to the standard inner product, 
such that 



H^=r?HrT 2 h = rjHr,- 1 = h) . 



(2) 



The virtue of such a conjugate pair h and H is that 
they possess an identical eigenvalue spectrum, because 
the Hamiltonians lie in the same similarity class. The 
reality of the spectrum is guaranteed, since one of the 
Hamiltonians involved, i.e. h, is Hermitian. The solu- 
tions of the corresponding time-independent Schrodinger 
equations H$ — e$ and h<f> = ecj) are then simply related 
as 



$ = 77 



(3) 



Let us first discuss the quantum brachistochrone problem 
for these types of systems. 



A. <f> — > 4> via u = $ — > $ via U 

Taking the initial state \ipi), the final state \ipf) to be 
orthonormal states of a Hermitian Hamiltonian system 
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2 



as well is the simplest situation to investigate. Equiv- 
alent would be to investigate the non-Hermitian sys- 
tem obtained by the similarity transformations ||2} and 
([5]) . Here we want to solve the quantum brachistochronc 
problem in a slightly less stringent form as studied in 
@, H, 0, S @|- Instead of solving ([1]), we just consider 
the physical relevant matrix element and seek the mini- 
mal time t =: t, such that a given transition probability 
is reached. This means for normalized initial and final 
states we solve the equation 



B. 



via U = $ — > $ via u 



!(</>, Kr,o)^; 



($/|t/(T,0)^ 



0, 



(4) 



for t with given constant < (3 < 1. Here u(t,t') and 
U(t,t') = ■q^ 1 u{t,t')ri are the time-evolution operators, 
which evolve a wavefunction from time t' to t associ- 
ated to the Hermitian and non-Hermitian system, respec- 
tively. In order to make the inner product involving the 
<!>s meaningful, we have to change the metric in the sec- 
ond expression in ((4]) as was argued in [ll] 



161 



or 



equivalently use a biorthonormal basis, see e.g. [17|, [lj 
Here we use the convention (^pf\Oipi) v :— (ipf\r] 2 04>i} for 
the 77-inner product with O being some operator. Clearly, 
for given time-evolution operators and final and initial 
states a real solution for r for all values of (3 does not 
always exist. Natural choices are for instance f3 = 1 or 
the maximum transition amplitude. 

Starting with time independent Hamiltonians, the 
problem is solved in a straightforward manner as we sim- 



ply have w(r,0) = e~ tTh and U(t,0) 



n -irH 



When 



cf>i and </>/ are orthogonal states we can find a solution 
in complete generality. Taking \<fr+) and \<f>-) to be two 
normalized eigenstates of h, the two states 



1 



|0 /A ) = ^(|0-}±*|0 + )). 



(5) 



are orthonormal with regard to the standard inner prod- 
uct. It is then straightforward to compute the matrix 
element occurring in ([4]) 



{<t>f\ 



-ith , 



1 



— ite- 



—ites. 



where the transition frequency between the two states 
is denoted as oj = e+ — e_. This means the quantum 
brachistochrone problem in the version Q is solved for 
the passage time 



t = — arcsin/3. 



(7) 



For P = 1 we recover r = tt/lj, which in slightly differ- 
ent forms, is well known and holds for any Hermitian or 
equivalent non-Hermitian system, as specified in As 
pointed out first in Q more spectacular results can be 
obtained when involving non-Hermitian Hamiltonians in 
the evolution operator while keeping the eigenstates to 
be associated to a Hermitian system. 



In [6(] a Gedankenexperiment was proposed in which a 
particle passes through a region, which causes its govern- 
ing Hamiltonian to change from a Hermitian to a non- 
Hermitian one. This scenario implies that the Hamilto- 
nian becomes explicitly time-dependent. The situation 
considered in [6J was for the initial and final states to 
be orthogonal states in a Hermitian system, whereas the 
time evolution was associated to a non-Hermitian Hamil- 
tonian. In general, we can write this temporary change 
of the Hamiltonian to a non-Hermitian Hamiltonian in 
the form 



H(t) =h + ghi(t), 



where h = and hi 7^ h\. This means we consider 
an analogue to a conventional time dependent scenario, 
however, with the difference that the perturbation is now 
non-Hermitian. A standard example of H{t) in ([5]) with 
hi being Hermitian is for instance the Stark-LoSurdo 
Hamiltonian describing an atom in an external electric 
field, with h representing the unperturbed atomic sys- 
tem and hi(t) the external electric field. In 0, Hfj] the 
alternative scenario was considered in which also the un- 
perturbed system was taken to be non-Hermitian. The 
treatment in [6| corresponds to the special case of ([H]) in 
which the time dependence is of the form of a stcpfunc- 
tion. This means to describe that setting one has to take 
H(t) = H = h + ghi for < t < t and H{t) = h for 
t > r, with t being the passage time rather than the 
pulse length as in the aforementioned example. 

Then, depending on the choice of the initial and final 
states, the time-evolution operator and the inner prod- 
uct, the quantum brachistochrone problem can be for- 
mulated in various different ways from (j?]). For instance, 
when projecting between orthogonal Hermitian states via 
a non-Hermitian time-evolution operator one may con- 
sider 



{<f, f \U(T,0)<lH 



\<t>f\\J<Pi 



I (^/Kr, 0)7^)1 
\v<t>f\ 



(g) where we used the ?7-norm defined as \\<f>\\. 



: P, (9) 

= Wl = 

\J {rj<t> I rj4>) . This switching between the eigenstates is in 
fact the key point of the entire analysis. We may view 
equation §§§ in two equivalent ways. On one hand on the 
right hand side we just project with a standard Hermitian 
time-evolution operator, between two somewhat unusual, 
but perfectly viable initial and final states rj4n and r]4>f 1 
respectively. This is a picture entirely in the standard 
quantum mechanical description with the difference that 
the initial and final states are not taken to be orthogonal. 
On the other hand we may use the equality and view this 
expression as a projection of some orthogonal initial and 
final states by means of a non-Hermitian time-evolution 
operator. We stress the fact that the initial and final 
states do not propagate in time and that the metric is 
not changed in this procedure [l9l |. 
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Slightly different and less natural is the possibility cor- 
responding to ([T]) when projecting with the standard in- 
ner product onto a final state. When written as an ex- 
pectation value this amounts to 



A 2 x 2 matrix Hamiltonian 



A pair of 2 x 2 matrix Hamiltonians related by a sim- 
ilarity transformation as in ^ is 



!((/>/ 1 c/(r, o)& 



(#/ \ u ( T i 0)#i 



= /8. (10) 



Alternatively, one may also envisage a situation when 
one projects from orthogonal Hermitian states onto 
eigenstates of a non-Hcrmitian Hamiltonian via a non- 
Hermitian time-evolution operator, or vice versa. Then 
one should solve 



[*f\U(T,0)<f>i 



$1 



| (0f\u(T, 0)^)| 

10/1 \ncj>i\ 



P, (ii) 



or 



{4> f \u{T,o)*i) v 1(^1^,0)^)1 



l„ll*4, 



/3- 



(12) 



As in a conventional time-dependent scenarios one is 
rarely able to compute the time-evolution operator ex- 
actly. However, assuming the non-Hermitian term in ([8]) 
to be small when compared with h, we may apply stan- 
dard pert urbation theory by iterating the DuHamel for- 
mula im, m, m 



U H {t,t') = U h {t,t')-ig / d8U H (t,a)hi(s)U h {a,l!). 



Taking Uh(t, t') = cxp(—ih(t—t')) we obtain to first order 
in g 



U H (t,0)=e 



— iht 



ig I ds e~ lh{t ' s) h 1 {s)e- lhs . (13) 
'o 



We may then compute for instance perturbatively the 
matrix element 



(<P f \U H (t,0)&) = (1-e-***) 



- +ig{(/> + \Imhi(j)-) 



When hi is Hermitian we naturally recover the result 
in ([6]). One may now proceed perturbatively using the 
above expression for Unit, 0). However, it is clear from 
the previous discussion that essentially all aspects of the 
problem, which we wish to consider here, may be illus- 
trated by selecting a two-level system from the larger, 
possibly even infinite, spectrum. Thus without loss of 
generality one may consider a 2 x 2 matrix Hamiltonian. 

In order to set the scene for the next section let us 
briefly recall with some minor variation the analysis of 



H 



re 
s 



Hi 



rcosQ 



2 

r cos 9 



(14) 



with lo = 2 V ' s 2 — r 2 sin 2 9 and r, s,9 € M. For the 
eigenvalues e± — rcos9 ± uj/2 to be real, one requires 
s 2 > r 2 sin 2 9, such that it is meaningful to introduce 
the new parameterization sin a = r/s sin 9 with a G M. 
This parameter range guarantees therefore unbroken VT- 
symmetry. The Hamiltonians h and H are related by the 
similarity transformation in @ , involving the Hermitian 
operator 



1 



sin a/2 —i cos a/2 
^/cosa \ i cos a/2 sin a/2 

The normalized eigenstates of H and h are 



(15) 



!$+)„ = 



e i 



^(1T1) 



01 V2 cos a V Te T 



I0d 



e ¥(l±i) ( ! 

Tl 



V2 



( 16 ) 

respectively. Taking now as initial and final states the 
orthogonal states \4>i) and as defined in (J5]), we com- 
pute 



ilk 



10* 



i6i 



COS f 

i sin f 



(17) 



and recover from J4|, with (3 — 1, the passage time 
t = tt/lo. This is what we expect from the general ex- 
pression ([7]) and in fact it is the same expression as ob- 
tained in 6|, where essentially the second equation in 
(j4|) was evaluated. On the other hand, if we now let 
the particle pass through the region in which the corre- 
sponding Hamiltonian becomes non-Hermitian, we may 
compute r by analyzing (jQj) or possibly (flQ|) . Acting with 



u(t,0) 



/III 



on the transformed initial state yields 



it) I 



»7|0i) 



itr cose / sin |( Q _ tuo) 



^/cos a \i cos | (a — tuj) 



(18) 



-ith 



When not acting on eigenstates with the operator e 
or e~ ltH one has to turn the infinite sum of operators 
into a matrix multiplication, see e.g. [f|. For this one 
can exploit the fact that any 2 x 2-matrix M can be 
decomposed in terms of Pauli matrices as M = ■ a 

with Hi € C, i = 0,1,2,3. Having expressed h or H 
in this manner, the operation with e~ uh or e~ ltH on a 
state reduces to a simple matrix multiplication by using 



the identity e ¥ 



cos c 



i smipfj, ■ a. 



Using (Tj"5]) . the matrix element in is computed to 



(<£/ 1^,0)04, = ^- 



i sin(a 



-) 



cos a 



(19) 
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Choosing the constant (3=1 and substituting ((19)) into 
©, we compute with \\<j)f\\ WfaWq = 1/ cosa the passage 
time t = ir/ui + 2a /ui. This expression involves now the 
parameter a, which may be tuned to make r arbitrarily 
small while keeping the transition frequency constant, as 
was first pointed out in [||. 

There are two equivalent ways of looking this result. 
On one hand we may think that one has solved the quan- 
tum brachistochrone problem entirely within the frame- 
work of the Hermitian system for some states, which have 
no obvious intrinsic meaning without referring to the 
non-Hermitian counterpart. On the other hand we may 
think that one has solved the time-dependent problem 
as outlined in the previous subsection involving the evo- 
lution with a non-Hermitian Hamiltonian between two 
orthonomal states in the Hermitian system. 

Similarly, we may consider a situation in which the 
final state is constructed from eigenstates of the non- 
Hermitian system and compute instead 



{$ f \U{T,0)ct> i ) v =ie 



i cos I (a — tlo) 



(20) 



^/cos a 

This yields for the same choice of fi = 1 the passage 
time r = 2n/uj + a/uj with ||$/|| H^ill,, = 1/V C0SQ! - We 
summarize our findings for the different types of scenarios 
in the following table: 



Hi) 




U(t,0) 


T 


1*0 




e - lH r 


tt/uj 


\<t>i) 


\h) 


e -iHr 


tt/lo + 2u/lo 


\<t>i) 


1*/} 


£ - lH r 


2tt/lo + a/u 






e -iHr 


2tt/lo + a/u 



Table 1: Passage times for various different initial and final 
states evolved by means of a PT-symmetric non-Hermitian 
time evolution operator. All matrix elements are computed 
with the 77-inner product and states are normalized with the 
77- norm. 

Next we demonstrate that such type of behaviour is 
not limited to PT-symmetric Hamiltonians, but can also 
be found for dissipative systems, i.e. genuinely non- 
Hermitian Hamiltonians with complex eigenvalues with 
negative imaginary part. 



III. NON-HERMITIAN DISSIPATIVE 
HAMILTONIAN SYSTEMS 

As argued above it is sufficient to consider a 2 x 2 
matrix Hamiltonian. We will now consider two different 
types of dissipative systems, i.e. those which have real 
and those with complex transition frequencies. 



A. Real transition frequency 

Let us modify the Hamiltonian H in (fT4")l slightly, such 
that it becomes a genuinely dissipative system. In order 



to achieve this we need to break the PT-symmetry not 
only for the wavefunction, but also for the Hamiltonian. 
Such type of Hamiltonians result for instance as effective 
Hamiltonians by coupling two non-degenerate states to 
some open channel as for instance explained in [24], [2f| . 
We consider here such type of Hamiltonian of the partic- 
ular form 



H 




iX 



s re 



-iO 



(21) 



with E, e, r, s,9,X € M. Similarly as for (jT4|) we will not 
provide here a concrete physical meaning for the parame- 
ters, as we would like to keep our treatment as generic as 
possible. Note that this Hamiltonian does not simply cor- 
respond to going to the regime of broken 'PT-symmetry 
for the Hamiltonian (fT4")l of the previous section. Instead, 
in the simultaneous limit E, e — ► and A — > i, the dissipa- 
tive Hamiltonian system H reduces to the PT-symmctric 
Hamiltonian H . 

The eigenvalues of H are computed to 



u 

E ± irX cos ( 

2 



(22) 



with u = e + — e_ = 2y/ (e + rX sin 9) 2 — X 2 s 2 denoting 
the transition frequency. For A being restricted to the 
interval — e/(s + sm9) < X < e/(s — sin 6), we can guar- 
antee that Co G M. In this parameter range the complex 
energy eigenvalues are indeed of the desired form of a 
decaying state, that is e± = E± — iT/2 with decay width 
T = 2rAcos6» <E M+ when -tt/2 < 6 < tt/2. It is then 
useful to introduce the parameterization 



sA 



sin a 



e + r X sin 6 ' 



(23) 



such that tana = 2sX/Q. The right eigenvectors of H 
corresponding to the eigenvalues e± in (|22p may then be 
expressed as 



1 



V2 



(l*->a±i !*+>*). 



(24) 



where |$±) a is defined in equation (fT6|) . From the consid- 
erations in the previous section it is clear that the opera- 
tor 77 is vital for the computations of the matrix elements 
occurring in the quantum brachistochrone problem, es- 
pecially when one wishes to evolve eigenstates of a Her- 
mitian Hamiltonian with a time-evolution operator asso- 
ciated to a non-Hermitian system. However, since for the 
case at hand the Hamiltonian H is now genuinely com- 
plex there can not exist any similarity transformation, 
which relates it to a Hermitian Hamiltonian. Nonethe- 
less, we can use the other property of 77, namely that it 
can be utilized to introduce a physically well defined in- 
ner product. This means we can seek a transformation 
such that the eigenstates fllMl) become orthonormal with 
regard to this product. From it is clear that we can 
take the same form for 77, but only have to replace a by 



5 



a to define a new fj. With the help of this new operator 
we construct the eigenstates \cf>±) — f]\$±), which yield 
indeed the desired orthogonality relations 



V 



(25) 

for n, m 6 {+,—}. The states 4> are eigenstates to 
the analogue of the Hermitian counterpart of a pseudo- 
Hermitian Hamiltonian. In fact, the adjoint action of f] 
diagonalizes H as 



h = fjHfj 1 = 











(26) 



Obviously we have now h ^ h\ but h has the same eigen- 
values as H, because it lies in the same similarity class. 



1. $ — > $ via U 

We are now in the position to solve the quantum 
brachistochrone problem for dissipative non-Hermitian 
Hamiltonians. For this we note first that when we try 
to compute the passage time r directly from the relation 
(HJ with (3 = 1, it will turn out to be complex. This 
would of course always be the case when the transition 
amplitude for all real values of t is smaller than f3. There- 
fore in order to find a physical solution we need to refor- 
mulate the quantum brachistochrone problem slightly to 
accommodate also the dissipative non-Hermitian Hamil- 
tonians. A natural expression to consider is one which 
features explicitly the decay width T, such as 



<t>f 



4> 



(27) 



for normalized initial and final states with < f3 < 1. 
This means for stable particles, i.e. r — > 0, we recover the 
expression . In analogy to the normalized eigenstates 
of the Hermitian Hamiltonian h in ([5]), we take now the 
initial and final states to be 



4>f/i 



V ( 
V2 V 



$_ > ±i 



-)±i 



Assuming at first no time dependence in the Hamilto- 
nian, we compute 



fa 



-ith ' 



l 



(l - e- lM ) . (29) 



Solving then ([27| for the passage time in complete anal- 
ogy to ([4]) with (3 — 1 yields the same expression for the 
passage time, namely r = tt/uj. Alternatively we could 
have also computed the second expression in ([27]) using 
U(t, 0) = e~ %m , which would of course lead to the same 
result. Possibly more interesting passage times can be 
obtained when we evolve the states (|2"5)) with the ana- 
logue to the non-Hermitian time evolution. 



2. $ — > <& via u 

We may assume now an explicit time dependence as in 
section II C and try to evolve the states by means 
of a time-evolution operator involving the Hamiltonians 
H in (|2ip. For that situation to make sense we need to 
consider the 77-inner products and generalize ^ to 



= f3e- 



(30) 

by introducing the decay width on the right hand side, 
similarly as we extended © to (27|). This way the pas- 
sage time will result to be real. Assuming a simple step- 
function time-dependence in ([8]), we compute 



{fa 


u(T,o)4>i) 




{vfa 


u(T,0)fj4>i) 




fa 


n 


4>i 


V 


vfa 


v4>i 





e- ith fj 



-it(E~i\r cos 9) ^/ cos 5, /_ e 'itcD/2 



V2(cosf +sinf ) 
and the matrix element 



-itui/2 



(j> f 



-it(E 



-iXrcosd) 



(31) 
(32) 

(33) 



Since = 1 in (|30p . we obtain with (3 — 1 the 

passage time f = ir/Q. Likewise we may evaluate ana- 
logues to the other matrix elements computed in section 
II C. Our results are summarized in table 2 



1M 


V'/) 


U(t,0) 


T 


4 


h) 


e- u " 


tt/u> 




fa) 








h) 


e~ u " 


tt/lj 




fa) 


e -itH 


tt/lj 



Table 2: Passage times for various different initial and final 
states evolved by means of a dissipative non-Hermitian time 
evolution operator with real transition frequency. All matrix 
elements are computed with the rji-inner product and states 
are normalized with the r;-norm. 

Thus somewhat surprisingly, despite the fact the ma- 
trix elements are somwhat different, the normalization 
factors compensate for this and in all cases the passage 
time results to f = n/Q. This might suggest that we 
should really attribute the possibility of tunable passage 
times to the 'PT-symmetry. However, we have not yet 
studied the possibility when uj ^ R. 



B. Complex transition frequency 

Often it is not even possible to restrict the parameter 
range so nicely like in the previous subsection as to ensure 



6 



that Cj £ R. Instead considering the scenario of leaving 
this regime for the above example, let us consider a dif- 
ferent system, which does not evenpossess such a regime 
and can be found for instance in 126] 



H 




2 cos </> sin 2(f) \ 
sin 2(f> 2 sin 2 <fi J ' 

_ (34) 
with E, e £ R and A, tj> £ C. In [26j the special case 
A, <p G R was treated for a concrete physical setting. Here 
we keep these parameters completely generic. The eigen- 
values of H are 



E± 



UJ 



A 



(35) 



the choice (3 > (3' would lead to unphysical complex pas- 
sage times. In order to keep matters simple, we could 
make a natural choice and just take this maximum of the 
right hand side of (f3"5| . i.e. (3 — (3' . However, due to 
the transcendental nature of this equation, there is no 
elegant analytic solution for r to this equation. Instead 
we try to make this solution to be as closely related to 
previously computed expressions as possible and choose 
the right hand side such that the passage time becomes 
t = n/Cj r leading to a particular p < (3'. We note that 
this choice does not lead to a loss of generality with re- 
gard to the main aim of our investigation, which is to seek 
passage times which can be made arbitrarily small. Tak- 
ing any other value below the maximum will simply lead 
to another definite value for r, when keeping Cj constant. 



with energy gap Cj — i+ — e_ = -\/4e 2 — A 2 — 4ieAcos2</>. 
Now we have lost the property of the transition frequency 
to be real. We parameterize instead 



-iXsin2(f> 



2e + Cj — i\ cos 20 



(36) 



such that the eigenvectors can simply be taken to be 
l^±) = |3?±)&- Replacing now also in 77, as defined in 
(j!5[) . the parameter a by a defines a new operator fj. We 
employ this operator to construct the Hamiltonian 



2. i-t$M« 

As in all previous scenarios we take the special case 
$ — > $ via U as a benchmark for the choice of the con- 
stant (3 = (3, since it can be dealt with analytically. Let 
us now compute 



(fa 






1 fj 



cosh(tCji) — cos(2d r — td r ) 
2 e Kt cos d cos d* 



(39) 



h = fjHfj- 1 = 



(37) 



with eigenvectors \<j)±) := \<j>±) as defined in (jTHJ) . Wc 
may now compute the passage time in a very similar fash- 
ion as for the 'PT-symmetric example, keeping however 
in mind that the transition frequency Cj as well as the 
parameter a are complex. A consequence of the latter is 
that r) is no longer Hermitian, that is ff ^ fj. For con- 
venience we introduce the abbreviations Cj = d r + iu>i, 
a = a r + ia.i and A = A,. + i\ for Aj/ r , Wj/ r , dj/ r £ R. 



With 



fa 



cosh a* 



cos a cos a 



(40) 



and the choice of (3 — $ as discussed in the previous 
section, namely taking it to be the right hand side of (|38|) 
with t — > -k/Cjj.^ the quantum brachistochrone problem 
amounts to solving 



cosh(ttDi) — cos(2d r — tCj r ) 
cosh 2 &i (1 + cosh(7r Cj j / ti r )) 



(41) 



$ — > $ via U 



It is straightforward to compute the square of the tran- 
sition probability 



fa 



-ith ' 



— e 
2 



-Ar* 



cosh(ia)j) — cos(tov)] , (38) 



when taking the initial and final states to be the or- 
thonormal states as defined in (JS)), with \<f>±) — > \4>±). 
There are now various possibilities we can equate this to 
and subsequently compute a passage time r. However, 
choosing a generic < (3 < 1 as for instance in equation 
(|4]) would lead to very involved solutions, due to the tran- 
scendental nature of this equation. Furthermore, since 
we are now dealing with a dissipative system the right 
hand side of {551) only reaches a maximum (3' < 1. Thus 



in this case. Since this is a transcendental equation, we 
can not solve it in complete generality and we are there- 
fore content to discuss some numerical solutions. For 
this purpose we can solve the equation of the transi- 
tion frequency for the coupling constant A = — 2i£cos<^+ 
\J 4e 2 sin 2 2<j) — to 2 and express it as a function of tfi, e and 
w. Since we want to keep u> constant, we investigate (|41[) 
as a function of time t by varying c\> and A or s and A see 
figure 1 or figure 2, respectively. 

The analysis of equation (|4Tj) . as depicted in figures 
1 and 2 demonstrates that it is possible to find tunable 
passage times for Hamiltonian systems of the type (|34|) . 
The precise dependence of the system on the parameters 
is rather involved in this case, but our analysis demon- 
strates that it is possible to approach r ~ 0. Similar 
conclusions can be drawn when changing \<f>f) to \$f) or 
\<f>i) to !<&,-), respectively. 
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0.25 
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-0.25 



0.2 0.3 

time 

Figure 1: Right hand side(rhs) minus the left hand side(lhs) 
of equation (141[) as a function of time t with fixed values of 
£ = 2.5, £ = 2,w = 4.87 - i3.31 for fa = 0.2, Ai = 3.5 ; 
02 = 0.3, A 2 = 3.73 + 20.2 ; <fe = 0.35, A 3 = 3.87 -M0.34 ; 
4>4 = 0.4, A4 = 4.02 + i0.52. The corresponding passage times 
t are simply taken to be the smallest values of t for which the 
rhs-lhs=0. 
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Figure 2: Right hand side minus the left hand side of 
equation (|41l) as a function of time t with fixed values of 
4>= 0.2, E = 2,w = 6.99-i0.46 for £1 = 9.0, Ai = 1.85-il4.84 
; £2 = 9.5, A 2 = 2.72 - il6.32 ; £3 = 9.9, A 3 = 3.41 - il7.29 ; 
£4 = 11.0, A 4 = 5.01 - il9.62. 



IV. CONCLUSION 

In Q the authors have extended the formulation of 
the quantum brachistochrone problem by allowing that 
the time-evolution operator may be associated to non- 
Hermitian Hamiltonians which are PT-symmetric, that 
is those with real eigenvalues. Here we did not insist on 
the PT-symmetry of the Hamiltonian, but allowed this 
symmetry to be completely broken. In order to take the 
most extreme case, we did not just spontaneously break 
the PT-symmetry for the wavefunction, which would re- 
sult in complex conjugate pairs for the energy eigenfunc- 
tion, but we allowed in addition that the PT-symmetry 
is also broken for the Hamiltonian. Thus we have consid- 
ered an effective Hamiltonian whose energy eigenvalues 
have a negative imaginary part, such that it is associated 
to dissipative systems. We found the same intriguing fea- 
ture as observed in Q for the quantum brachistochrone 
problem for PT-symmetric non-Hermitian Hamiltonians, 
namely that the passage time can be made arbitrarily 
small also for non-Hermitian Hamiltonians associated to 
these type of systems. Our observations suggest that 
this type of phenomenon may occur when one projects 
between orthonormal states, which are not eigenstates of 
the non-Hermitian Hamiltonian associated to the time- 
evolution operator, irrespective of whether this Hamilto- 
nian is PT-symmetric or not. 

Clearly there are various open questions to be an- 
swered. First of all, it would be highly desirable to 
have a more formal and generic proof for this phenom- 
ena, rather than case-by-case studies. This holds for the 
PT-symmetric case treated in Q as well as for the case 
presented here. In addition, one could make the above 
considerations more involved by allowing more compli- 
cated time-dependences rather than the simple stepfunc- 
tion and study other possibilities in (J8J) . For such more 
realistic scenarios we may have to resort to a perturbative 
treatment using (flT 
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